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ABSTRACT 

CD \ We had reported in IChatv et al.l l)2002a(l on the near-infrared (NIR) identification of 

^vO ■ a possible counterpart to the black hole candidate XTE J1908+094 obtained with the 

' ESO/NTT. Here, we present new, follow-up, CFHT adaptive optics observations of 

the XTE J1908+094 field, which resolved the previously proposed counterpart in two 
objects separated by about 0"8. Assuming that both objects are potential candidate 
counterparts, we derive that the binary system is a low-mass system with a companion 
star which could be either an intermediate/late type (A-K) main sequence star at a 
' distance of 3-10 kpc, or a late-type (>K) main sequence star at a distance of 1-3 

kpc. However, we show that the brighter of the two objects (J ~ 20.1, H ~ 18.7, 
q . K' ~ 17.8) is more likely to be the real counterpart of the X-ray source. Its position 

!h ' is more compatible with our astrometric solution, and colours and magnitudes of the 

, other object are not consistent with the lower limit of 3 kpc derived independently from 

d • the peak bolometric flux of XTE J1908+094. Further multi-wavelength observations 

of both candidate counterparts are crucial in order to solve the pending identification. 
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1 INTRODUCTION 

The X-ray transient XTE J1908+094 was serendipitously 
discovered on 2002 February 21 during observations of the 
field of the Soft Gamma-ray Repeater SGR 1900+14 per- 
formed with the Proportional Counter Array (PCA) aboard 
RXTE (galactic c oordinates l,b = 43.26 deg, +0.43 deg; 
IWoods et all 120021 Within one month after its discov- 
ery, the source flux increased by a factor ~ 3 and broad 
peak QPOs at a frequency of 1 Hz were also detected. 
The RXTE/PCA 2-30 keV spectrum was found consis- 
tent with a power-law (photon index=1.55) w ith an hydro- 
gen c olumn density Nh = 2.3 x 10 22 cm -2 iWoods et alJ 
2002). Soon after, XTE J1908+094 was observed at higher 
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energies (15 — 300 ke V) with the PDS aboard BeppoSax 
jFeroci fe Reboall2002lh at a flux level inconsistent with the 
extrapolation of the RXTE/PCA power law, implying a 
cut-off in the X-ray spectrum at energies around 100 keV. 
Target of Opportunity ( ToO) observations perfor med with 
the MECS on BeppoSax Jin t Zand et al.ll2002ah provided 
a source positioning better than 20", while the measure- 
ment of the source flux in outburst set a lower limit 
on the distance of 3 kpc Jin't Zand et alJ l2002blh The 
transient radio counterpart to XTE J1908+094 was iden- 
tified with the Very Large Array (VLA), with a posi- 
tion uncertainty of 0.1" in each coordinate iRupen et alJ 
|2002j) , consistent with t he improved Chandra/ ACIS position 
Jin't Zand et all2002p|) . First follow-up optical observations 
fail ed to detect any possible c ounterpart brighte r than R ~ 
23 iWaener fc Starrfieldl 2002) and I ~ 22 iGarnavich et all 
120021) . Thereafter. IChatv et alJ (1200211 reported the discov- 
ery of a likely near-infrared (NIR) counterpart, appare ntly 
variable along the X-ray lightcurve llChatv et alJl2002aD . 
Here, we report on the results of new, follow-up, adaptive 
optics NIR observations of the XTE J1908+094 candidate 
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counterpart. The observations are described in Section [5] 
while the results are discussed in Section [3] 



2 OBSERVATIONS AND DATA REDUCTION 

The first NIR observations of the XTE J1908+094 field were 
obtained in April and June 2002 with the SOFI camera at 
the ESO 3.5 m New Tech n ology T elescope (NTT) and are 
described in IChatv et al.l J2002a). New NIR observations 
were performed on 2002 August 18 in visitor mode with 
the 3.6m Canada-France-Hawaii Telescope (CHFT) on top 
of Mauna Kea (Hawaii). The observations were carried out 
using the Adaptive Optics Bonnette (PUEO), mounted at 
the f/8 CFHT Cassegrain focus, equipped with the KIR in- 
strument. This is a near-infrared imaging camera based on 
the 1024x1024 Rockwell Science Center HAWAII (HgCdTe 
Astronomical Wide Area Infrared Imaging) focal plane ar- 
ray, sensitive to radiation from 0.7 to 2.5 /im. The camera 
plate scale is 0"035/pixel, yielding a total field-of-view of 
36" x 36". 

Images were taken through the J (A = 1.250 /im, 
FWHM = 0.160 fan), H (A = 1.635 fan, FWHM = 
0.290 fan) and K' (A = 2.120pm, FWHM = 0.340pm) fil- 
ters. The source was observed through a sequence of 5 ex- 
posures of 90s each in the J, H and K' filters, respectively, 
with a random offset among images in order to perform back- 
ground subtraction of the variable IR sky. Conditions were 
photometric. The airmass was between 1.02 and 1.04, the 
typical seeing conditions were around 0.4" — 0.6" and the 
adaptive optics resulted in a source Point Spread Function 
(PSF) of about 0"1-0"14 in the images (diffraction limited). 
The images have been corrected for instrumental effects, by 
removal of the dark current and flat-fielding, and then sky 
subtracted, using standard procedures available in the IRAF 
package ccdred. For each filter, we finally computed a final 
image which is the median of the sequence. Photometric cal- 
ibration was perform ed using the standard star field FS35 
fG21-15: ICasalll992Tl . 



3 RESULTS 
3.1 Astrometry 

To accurately regis ter the radio coordinates of the source 
jRupen et al] £20021 on the CFHT images, we have com- 
puted the image astrometry using as a reference the posi- 
tions of stars selected from the Two Micron All Sky Survey 
(2MASS) Point Source Catalogue, which has an intrinsic ab- 
solute astrometric accuracy of ~ 0"2 per coordinate 1 . Due 
to the smaller field of view of the CFHT camera (36" x 36"), 
only 11 2MASS objects have been identified in the averaged 
J-band image and used as suitable astrometric calibrators. 
The pixel coordinates of the reference stars have been com- 
puted by a two-dimensional gaussian fitting procedure, and 
transformation from pixel to sky coordinates was then com- 
puted using the programme ASTROM 2 , yielding an rms of 
~ 0"08 in both Right Ascension and Declination, which we 

1 http: / /spider. ipac.caltech.edu/staff/hlm/2mass/overv/overv. html 

2 http:/ /star-www. rl.ac.uk/Software/software. htm 



assume representative of the accuracy of our astrometric so- 
lution. After taking into account the rms of the astrometric 
fit (~ 0"08), the absolute accuracy of the 2MASS coordi- 
nates (~_0j^2j_jmd th^error of the source radio coordinates 
fO"l. lRupen et aklEofB) . we estimated ~ 0"24 (per coor- 
dinate) to be the final uncertainty on the computed source 
position. 

Figure □ shows 15" x 15" and 4"5 x 4"5 cutouts of the 
NTT and CFHT J-band images centered around the com- 
puted source position (left and right panel, respectively). In 
the CFHT image we overplot the 1 and 3 a error circles de- 
rived from our astrometric solution. One and two objects, 
labelled 1 and 2 in figure (right panel), are found within 
the la and 3a error circles, respectively. As the angular sep- 
aration between objects 1 and 2 is about 0"8, we note that 
they could have bee n hardly resolved in the SOFI images of 
IChatv et aI1j2002al) . in view of the non-optimal seeing con- 
ditions (0"8 — l"2) during the observations and of the much 
larger pixel size of the detector (0"292). Thus, we conclude 
that the previously proposed candidate counterpart was ac- 
tually the blend of these two objects. 

Although both objects are considered as potential coun- 
terparts in the following analysis, we note that object 1 (the 
brighter of the two) is certainly favoured as its position is 
more consistent with our astrometric solution. 

3.2 Photometry 

We have computed the JHK' photometry of the two candi- 
date counterparts (objects 1 and 2) identified in the CFHT 
images. In order to obtain more precise measures, we have 
used the IRAF/DAOPHOT package for the photometry in 
crowded fields. Specifically, we computed a typical PSF for 
every image using the tasks daofind, phot, psf, nstar and sub- 
star, iterating the procedure to minimize the residuals after 
PSF subtraction. Then, we ran a star subtractio n algorithm 
through the task allstar iMassev fc Davislll992f) to remove 
all objects detected in the neighbourhood of our candidate 
counterparts before computing their magnitudes. The final 
CFHT JHK' magnitudes of objects 1 and 2 are reported in 
rows 4 and 5 of Table 

Using the same procedure describ ed above, we have re- 
analyzed the original NTT images of IChatv et al] (l2002al) 
in an attempt to resolve, or at least constrain, the flux of 
objects 1 and 2. Unfortunately, the highly uncertain resid- 
uals after the PSF subtraction make it impossible to de- 
rive reliable magnitude estimates for b oth objects. Thus, w e 
have taken the magnitudes reported in IChatv et ail (l2002af) . 
which correspond to the blend of objects 1 and 2, as upper 
limits on the actual magnitudes of these two objects. These 
values are listed in the first three rows of Table 

As no direct evaluation of the variability of objects 1 
and 2 is possible from the current photometry, we have tried 
to obtain an indirect piece of evidence by simply compar- 
ing the combined flux of the two objects measured in the 
CFHT images (Table □ last row) with the flux of the blend 
measured in the NTT ones (Table Q first three rows). As it 
is seen from Tabled in the H and K bands the combined 
flux of objects 1 and 2 is significantly fainter than the flux 
of the blend at the different epochs (it is also true in the 
J-band, however we point out that the large error attached 
to the 2002-06-17 measurement hide any variation). In other 
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Figure 1. Left panel: 15" X 15" cutout of the NTT/SOFI J-band image of the XTE J1908+094 field iChatv et alj|2002al) . with the 
previously suggested candidate counterpart marked by an arrow. North to the top, East to the left. Right panel: 4"5 X 4"5 cutout 
(corresponding to the box in the left panel) of the combined CFHT/KIR J band image of the same field. In the CFHT image, the 
two circles (of radius 0"24 and 0"71, respectively) correspond to the lcr and 3cr uncertainties on the source position deri ved from our 
astrometric solution and the target radio position (a(J2000) = 19 h 08 m 53!77 ± 0f007; <5(J2000) = +09°23'04!'9 ± O'.'l) from lRupen et alJ 
J2002}) . The two possible candidate counterparts detected within the 3<r error circle are labelled (1 and 2). 

words, the IR flux of one of the two objects (or both) varied 
in the time span covered by our observations. 

To search for a possible correlation between the ob- 
served IR and X-ray variability, we have compared our 
photometry reported in Table □ with the RXTE/ASM 
lightcurve of XTE J1908+094 (Figure H top panel). First, 
we considered as before the total JHK magnitudes of ob- 
jects 1 and 2 (Table 1, last row) and the magnitudes of 
the blend (Table 1, first three rows). These are plotted in 
Figure 2 (middle panel). Second, we considered the actual 
JHK magnitudes of objects 1 and 2 (Table 1, rows 4 and 5) 
with the expected magnitudes derived from the NTT/SOFI 
upper limits. As a crude approximation, we have straightly 
computed (per each passband) the difference between the 
NTT/SOFI flux of the blend and the CFHT/KIR flux of 
objects 1 and 2 in turn, assuming the other constant. The 
resulting lightcurves are shown in Figure [5] (lower panel). 
From the IR lightcurves it is clear that either object 1 or 2 
(or both) have undergone a significant decrease of the NIR 
flux consistent with the X-ray flux decay which occurred in 
the time span covered by our observations. This variability 
means that, almost certainly, one of the two objects is the 
real counterpart of the X-ray source XTE J1908+094. 

3.3 The companion star 

We have used our JHK' band photometry to derive con- 
straints on the nature of the two XTE J1908+094's candi- 
date counterparts. For both objects 1 and 2, we have taken 
the magnitudes measured on August 18th 2002, i.e. likely at 
the minimum, as inferred from the X-ray lightcurve (Figure 
[5] upper panel), assuming that all the flux is produced by 



the companion star. However, we have to caution that this 
assumption might not be entirely valid, since Figure [2] shows 
that the X-ray source still exhibited some level of activity 
at this epoch. Therefore, it is likely that the infrared emis- 
sion was not produced purely from the companion star but 
was contaminated by the contribution of the accretion disc. 
For this reason, we have prudently assumed the derived ab- 
solute magnitudes as upper limits on the luminosity of the 
companion star. 

We show in Figure [3] a theoretical J — K' , M K i, 
colour-magnitude diagr am (CMD) built for template stars 
fo.uelas-Mavorgal Il99ll) . We overplot the locations corre- 
sponding to the two XTE J1908+094 candidate counter- 
parts, with the corresponding absolute magnitudes being 
computed for ten different values of the distance, ranging 
from 1 to 10 kpc, and dereddened for three different val- 
ues of the hydrogen column densit y around the value of 
N H = 2.3 x 10 22 cm" 2 measured bv iWoods et alJ lj2002h . 

Since each derived value of the absolute magnitude rep- 
resents an upper limit, we can not accurately determine the 
distance of the source from the locations of their candidate 
counterparts on the CMD. However, the spectral type can 
be determined by the J — K' colour, which is independent of 
the source distance. Therefore, according to the diagram, the 
companion star of object 1 would be an intermediate/late 
type main sequence star, probably of spectral type between 
A and K, and at a distance between 3 and 10 kpc. In the 
case of object 2, the companion star would be a late-type 
main sequence star, probably of spectral type later than K, 
and at a distance between 1 and 3 kpc. 

The derived constraints on the spectral classification 
of the companion star are fully consistent with the non- 
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Date 


MJD 


Instrument 


Object 


J 


H 


K' 


2002.04.25 
2002.04.29 
2002.06.17 


52389.9 
52393.9 
52442.2 


NTT/SOFI 
NTT/SOFI 
NTT/SOFI 


Objects 1,2 
Objects 1,2 
Objects 1,2 


> 19.29 ±0.13 

> 19.86 ± 0.43 


> 17.21 ±0.21 

> 17.39 ±0.16 

> 17.86 ±0.07 


> 16.40 ±0.16 

> 16.52 ± 0.16 

> 16.65 ± 0.10 


2002.08.18 
2002.08.18 


52504.3 
52504.3 


CFHT/KIR 
CFHT/KIR 


Object 1 
Object 2 


20.14 ±0.14 
21.26 ±0.35 


18.69 ± 0.04 
18.98 ± 0.05 


17.77 ±0.03 
18.06 ±0.04 


2002.08.18 


52504.3 


CFHT/KIR 


Objects 1+2 


19.81 ±0.37 


18.07 ±0.07 


17.16 ±0.05 



Table 1. Photometry of the XTE J1908±094 candidate counterparts. The magn itudes in rows 1-3 w ere obtained with the SOFI 
instrument at the NTT and refer to the original candidate counterpart reported by IChatv et alj J2002al) . which is actually a blend of 
two objects clearly resolved by the CFHT (objects 1 and 2, see Figure^. Thus, the SOFI magnitudes have to be taken as upper limits 
on the magnitudes of both objects 1 and 2. The magnitudes in row 4-5 have been obtained with the KIR camera at the CFHT and refer 
to objects 1 and 2, respectively. Finally, the magnitudes in the last row represent the combined magnitudes of objects 1 and 2 (i.e. an 
artificial blending of CFHT objects) in order to compare it to the SOFI magnitudes. 




5.235xl0 4 5.240x10'' 5. 245*10'' 5.2S0xJ,0 4 5.255xl0 4 

MJD Time (JD - 2,400.000.5) 



Figure 2. Multi-wavelength lightcurve of XTE J1908±094. Upper panel: RXTE/ASM 2-10 keV lightcurve of XTE J1908±094. Flux 
values are in units of counts/s averaged over one day. Middle panel: observed IR lightcurves from the J (*), H (o) and K' (A) band 
magnitudes listed in Table Per each passband, the points at the first three epochs are the magnitudes of the blend of objects 1 and 
2 measured with the NTT/SOFI, while the points at the last epoch are the total magnitudes of objects 1 and 2 measured with the 
CFHT/KIR. Lower panel: Same as above, but the couple of points at the first three epochs are the estimated NTT/SOFI magnitudes for 
objects 1 and 2, corresponding to the difference between the total flux of the blend and the flux of object 2 and 1 (assumed constant), 
respectively. The couple of points at the last epoch are the measured CFHT/KIR magnitudes of objects 1 and 2. The dashed and dotted 
lines connect points associated to objects 1 and 2, respectively. The timescale is the same in the three panels, and the flux scale is the 
same in the middle and lower panels, for an easier comparison. 



det ection at R ^ 23 iWaener fc StarrfieldEool and at I ^ 
22 ijOarnavich et al]l2002f) . since a K or M type star at a 
distance ^ lkpc would yield, for the expected absorption, 
V ^ 30. Near- infrared observations are therefore of prime 
importance to reveal highly absorbed sources. 



However, taking into account the peak bolometric flux 
of XTE .11908+094. lin't Zand et alJ l)2002bl, have yielded a 
lower limit to the distance of 3 kpc (and even 6 kpc if the 
compact object is a black hole, with a minimum mass of 
3 Mq). This makes it unlikely that object 2 is the actual 
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Figure 3. J — K\ M K i , absolute CMD diagram computed for 
template main sequence stars iRuelas-Mavorgall99j) . The range 
of locations for the XTE J1908+094 companion star in the di- 
agram has been derived from the absolute magnitudes of the 
two possible counterparts, computed taking the minimum ap- 
parent magnitudes on 2002 August 18th, and assuming that 
all the flux is produced from the companion star. Three differ- 
ent values of the interstellar extinction have been applied i.e. 
Ay = 11.19,12.9,13.9 magnitudes from right to left, respec- 
tively. Ay = 12.9 m agnitudes corresp onds to a column density 
N H = 2.3 X 10 22 (Woods et alJl2002l) . using the relation from 
iPredehl fc Schmit j ll995l) . From bottom to top the crosses corre- 
spond to a source distance increasing from 1 to 10 kpc. According 
to the diagram, in the case of object 1, the companion star would 
be be an intermediate/late type main sequence star of spectral 
type between A and K, located between 3 and 10 kpc. In the 
case of object 2, the companion star would be a late-type main 
sequence star of spectral type later than K, located between 1 
and 3 kpc. 

counterpart of XTE J1908+094. Therefore, our analysis of 
the CMD favours, as does the astrometry, object 1 as the 
most likely counterpart. 



4 CONCLUSION 

We have reported on CFHT adaptive optics observations 
of the black hole candidate XTE J 1908+094 which showed 
that the previously reported counterpart llChatv et alJ 
2002a) was in fact a blend of two objects. A new and more 
accurate astrometric analysis provides qualitative evidence 
that the brighter of the two (called object 1) is a more likely 
counterpart. Accurate photometry performed on four sets 
of observations spanning 4 months have shown a decrease of 
the NIR flux for one (or both) of the two candidates, consis- 
tent with the X-ray decay, showing that one of the two is al- 
most certainly the real counterpart of XTE J 1908+094. Af- 
ter dereddening and assuming that the flux represents upper 
limits to the contribution of the companion star, we derived 
that the companion star could be i) either an intermedi- 
ate/late type main sequence star of spectral type between A 
and K, located between 3 and 10 kpc, or ii) a late-type main 
sequence star of spectral type later than K, located between 
1 and 3 kpc. However, since a lower limit to the distance 



of 3 kpc has been de rived from the peak bo lometric flux 
of XTE J 1908+094 bv lin't Zand et alJ J2002bl) . both colors 
and magnitudes, as well as the astrometry, would rule out 
object 2 and favour object 1 as the most likely counterpart 
of the X-ray source. 

Further near-infrared observations are needed, espe- 
cially during the next outburst of the source, in order to 
monitor the variability of both candidate counterparts and 
eventually confirm (or not) our identification of object 1 
as the counterpart of XTE J1908+094. Getting some spec- 
troscopy, although difficult because of faint and blended can- 
didates, should be feasible with next generation instruments, 
and would finally better constrain both the spectral type and 
distance of the companion star. 
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